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Abstract 

We study the use of the System Optimum (SO) Dynamic Traffic As¬ 
signment (DTA) problem to design optimal traffic flow controls for free¬ 
way networks as modeled by the Cell Transmission Model, using variable 
speed limit, ramp metering, and routing. We consider two optimal control 
problems: the DTA problem, where turning ratios are part of the control 
inputs, and the Freeway Network Control (FNC), where turning ratios 
are instead assigned exogenous parameters. It is known that relaxation 
of the supply and demand constraints in the cell-based formulations of 
the DTA problem results in a linear program. However, solutions to the 
relaxed problem can be infeasible with respect to traffic dynamics. Pre¬ 
vious work has shown that such solutions can be made feasible by proper 
choice of ramp metering and variable speed limit control for specific traffic 
networks. We extend this procedure to arbitrary networks and provide 
insight into the structure and robustness of the proposed optimal con¬ 
trollers. For a network consisting only of ordinary, merge, and diverge 
junctions, where the cells have linear demand functions and affine supply 
functions with identical slopes, and the cost is the total traffic volume, 
we show, using the maximum principle, that variable speed limits are not 
needed in order to achieve optimality in the FNC problem, and ramp me¬ 
tering is sufficient. We also prove bounds on perturbation of the controlled 
system trajectory in terms of perturbations in initial traffic volume and 
exogenous inflows. These bounds, which leverage monotonicity proper¬ 
ties of the controlled trajectory, are shown to be in close agreement with 
numerical simulation results. 
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1 Introduction 


The System Optimum (SO) Dynamic Traffic Assignment (DTA) problem, in¬ 
troduced in [201EU, has attracted significant interest from the transportation 
research community, see, e.g., |23| for an overview. While originally proposed 
mainly for planning purposes, it is also being increasingly used as a framework to 
compute optimal control for traffic flow over freeway networks, e.g., see mm, 
when traffic controllers aim to minimize a global cost of the whole network - 
hence Social Optimality, as opposed to the single-vehicle oriented User Equilib¬ 
rium modeling frameworks. Continuing along these relatively recent trends, this 
paper focuses on the use of solutions of two variants of the SO-DTA to design 
optimal controls for dynamic network traffic flows over a given time horizon, in 
the form of variable speed limit, ramp metering, and routing (turning ratios) 
matrices. 

The Cell Transmission Model (CTM), originally proposed in |8l[9], is a com¬ 
pelling framework to simulate realistic first order traffic dynamics. It consists 
of a time and space discretization of the kinematic wave models of Lighthill- 
Whitham and Richards [TSl El] ■ Unfortunately, straightforward formulations of 
DTA for the CTM are known to lead to non-convex problems, and hence are 
unsuitable especially for real-time applications. On the other hand, in the DTA 
formulation of m, the supply and demand constraints of the CTM are relaxed 
to yield a linear program. However, the computational simplicity resulting from 
this relaxation comes at the expense of possible infeasibility of a resulting op¬ 
timal solution with respect to traffic dynamics. Quite interestingly, |22j shows 
that the optimal solution of a linear program analogous to the DTA relaxation 
in [27] can be realized exactly for traffic dynamics modeled by the link-node cell 
transmission model by proper design of ramp metering and variable speed limit 
controller, when demand functions are linear, supply functions are affine, and 
the network consists of a mainline with on- and off-ramps. In this paper, we 
consider extensions of the approach proposed in mm to arbitrary networks, 
where traffic dynamics is inherited by the CTM, with traffic dynamics encom¬ 
passing the ones originally proposed in and also allowing for arbitrary 

concave demand and supply functions, and convex cost functions, including to¬ 
tal travel time, total travel distance, and total delay as special cases. Under 
these generalizations, the resulting SO-DTA is a convex program, which can be 
solved using readily available software tools such as cvx [T] [13] and possibly 
suited for distributed iterative solvers, e.g., see our preliminary work [1]. 

In short, the considered approach to the optimal control of freeway networks 
over a finite time horizon consists of two steps: (i) to formulate and solve convex 
optimal control problems that are relaxations either of the DTA problem (with 
turning ratios part of the optimization), or of the Freeway Network Control 
(FNC) problem (where the turning ratios are exogenously imposed); (ii) to 
design open-loop variable speed limits, ramp meters and routing controls over 
the time horizon to make such optimal solution feasible with respect to traf w 
dynamics modeled by the CTM. Natural questions concern: (a) under which 
conditions step (ii) above is not necessary, i.e., the optimal solution of the convex 
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optimal control problem is readily feasible with respect to traffic flow dynamics; 
and (b) how robust the optimal control computed through the procedure above 
is with respect to perturbations to the network during the execution of the open- 
loop controller over the time horizon. The main novel contributions of this paper 
address questions (a) and (b) as follows. On the one hand, using Pontryagin’s 
maximum principle, we prove that, for networks consisting only of ordinary, 
merge, and diverge junctions, and whose cells have linear demand and affine 
supply functions with identical slopes, the optimal solution of the FNC problem 
with total traffic volume as cost does not require the use of variable speed limits 
as a proper choice of ramp metering controls makes it readily feasible with 
respect to the CTM model of [5] . On the other hand, we provide bounds on the 
perturbation to the system trajectory under the open-loop controller obtained 
from the two-step procedure due to perturbations in the initial traffic volumes 
and exogenous inflows. In order to derive such bounds, we use the fact the 
nominal (i.e., unperturbed) controlled system trajectory resulting from steps (i) 
and (ii) above is always in free-flow, and hence it satisfies a certain monotonicity 
property that in turn can be used to evaluate its robustness to perturbations of 
the initial traffic volumes and exogenous inflows. 

It is helpful to clarify the novelty of our contributions with respect to existing 
literature. We generalize the applicability of the two-step procedure of using 
solution of the two DTA variants for design of optimal traffic flow control to 
general network topologies, concave supply and demand functions, convex cost 
functions, than the ones considered previously, e.g., in [121122j . The maximum 
principle has been used to identify necessary conditions for optimal control of 
traffic flow over networks, e.g., in m- However, the underlying model for traffic 
flow dynamics in m does not capture backward propagation of congestion, and 
in particular does not resemble CTM. 

Robustness of open-loop controllers can be quantified through standard sen¬ 
sitivity analysis of controlled traffic dynamics. Our bounds, which exploit the 
monotonicity properties of controlled system trajectories, are applicable to rel¬ 
atively larger perturbations than those obtained through such standard tech¬ 
niques. Implications of such monotonicity property for robustness of dynamic 
network flows have been recently investigated in different contexts [1 0 E]. 
Our robustness analysis of the solution to deterministic DTA problems is to be 
contrasted with chance-constrained solution of stochastic SO-DTA, e.g., in [26] 
under probabilistic information about inflows. Our approach to applying max¬ 
imum principle and robustness bounds necessitates consideration of continuous 
time versions of CTM in the analytical part of the paper. Such continuous time 
and discrete space versions have been used previously, e.g., to develop proba¬ 
bilistic versions of CTM [16] . Our adoption of continuous time version is merely 
to facilitate analysis, and is not to be interpreted as a new numerical framework 
for traffic flow dynamics. Indeed, the simulations reported in this paper are 
performed in the standard discrete time discrete space version of CTM. 

The remained of the paper is organized as follows. In Section the DTA 
and FNC are formulated as optimal control problems for continuous-time cell- 
based dynamic traffic models. In Section]^ convex optimal control relaxations 
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Figure 1: A multi-source multi-sink transportation network topology. 

of the DTA and FNC problems are presented and proved to be tight through 
a proper choice of open-loop variable speed, ramp metering, and routing con¬ 
trollers (Proposition . In Section scenarios under which no variable speed 
and ramp metering control is required in order to achieve optimality are identi¬ 
fied (Proposition]^. Sectionpresents perturbations bounds in optimal DTA 
and FNC solutions due to perturbations in initial traffic volumes and exogenous 
inflows: the main results are stated in Propositions andSection [O] reports 
some numerical simulations illustrating the main theoretical results. Finally, 
Section]^ draws conclusions and suggests future research directions. The proofs 
of the technical results are presented in the Appendix. 

2 Dynamic Traffic Model and Problem Formu¬ 
lation 

We describe the topology of the transportation network as a directed multi¬ 
graph Q = {V,£) with links i G £ representing cells and nodes representing 
either junctions (being them of merge, diverge, or mixed type) or interfaces 
between consecutive cells (briefly referred to as ordinary junctions). The head 
and tail nodes of a cell i are denoted by ti and at, respectively, so that the cell 
i = is directed from to r^. One particular node w GV represents the 

external world, with cells i such that ai = w representing sources (identifiable 
with on-ramps in freeway networks) and cells i such that ti = w representing 
sinks (representing off-ramps in freeway networks). The sets of sources and 
sinks are denoted by TZ and S, respectively. The network topology is typically 
illustrated by omitting such external node w and letting sources have no tail 
node and sinks have no head node. (See Figure]^) The notation 

= {(b j) G £ X £ : Ti = aj ^ w} 

is used to refer to the set of all pairs of adjacent (consecutive) cellsQ 

^Each cell is meant to model a portion of an actual road of a given length, so that the 
graph Q need not in general coincide with the actual road network topology. In fact Q may 
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Figure 2: Concave (left) and piecewise linear (right) fundamental diagrams. 
Observe that the set of pairs {x,z) such that 0 < z < min{d(a:), C, s(a;)} is a 
closed convex set (in the special case in the right figure it is in fact a polytope). 


The dynamic state of the network is described by a time-varying vector 
x{t) € whose entries Xi{t) represent the traffic volume in the cells i G £ at 
time t > 0. The inputs to the network are the exogenous inflows Xi{t) > 0 at the 
sources i G TZ, while the outputs are the outflows /ii(t) > 0 at the sinks i G S. 
Conventionally, we set Xi{t) = 0 for all non-source cells i G S\TZ and /ii(t) = 0 
for all non-sink cells i G £\S, and stack up all the inflows in a vector X{t) G 
and all the outflows in a vector /i(t) G We denote the traffic flow from a 
cell i to an adjacent cell j by the total inflow in and the total outflow 

from a cell i by yi{t) and Zi{t), respectively, and stack up all cell-to-cell flows, 
total inflows, and total outflows in vectors f{t) G y(t) G and z{t) G K^, 
respectively. The law of mass conservation then reads 

x^{t) = yi{t) - Zi{t), iG£, (1) 

yi{t) = X^{t), Zi{t) = n^{t) + , iG£, (2) 

je£ 

In fact, equation 0 states that the time-derivative of the traffic volume Xi 
on a cell equals the imbalance between its inflow yi and its outflow Zi, while 
equation ([^ states that the inflow yi in a cell is the aggregate of the exogenous 
inflow Xi and the flows fji from other cells in the network, while, symmetrically, 
the outflow Zi from a cell is the aggregate of the flows fij to other cells in the 
network and the outflow fj,i towards the external world0 

represent a refinement of such actual road network topology, with a given road possibly split 
into multiple adjacent cells. Observe that such refinement maintains the original junction 
layout as well as the original sets of sources and sinks. 

^In fact, one could have replaced yi and Zi in both the right-hand sides of 0 and 
by the expressions in the right-hand sides in 0 and reduced the number of variables in 
the DTA problem. However, introducing the variable yi and Zi turns out to be useful for the 
problems discussed in Section]^ Using these additional variables also proves useful in deriving 
distributed optimization algorithms based, e.g., on the alternating method of multipliers or 
interior point methods, e.g., see U. 
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Figure 3: Demand control is actuated by ramp metering on the sources i G TZ 
(left figure) and variable speed limit in the other cells i G £\TZ (right figure). 
As proved in Proposition]^ for any feasible pair {x*, z*) such that 0 < z* < 
min{(i(a;*), (7, s(a:*)} there exists a choice of the demand control parameter a 
such that z* = d{a,x*). 


We consider first order dynamic traffic models and assume a concave Fun¬ 
damental Diagram, as the one shown in Figure]^ to model the relationship be¬ 
tween traffic volume and flows. Following Daganzo’s Cell Transmission Model, 
we introduce demand and supply functions, which return the maximum outflow 
from and, respectively, the maximum inflow in a cell as a function of its current 
traffic volume and, possibly, additional control parameters. In turn, they can be 
interpreted as representing the rising and, respectively, descending parts of the 
Fundamental Diagram (see again Figure]^. In particular, the demand functions 
are assumed to take the form 

dt{xi,ai) = inm{aidi{xi),Ci}, iG£\TZ, (3) 

di{xi,at) = miii{d^{xi),aiQ} , iGlZ, (4) 

where: di{x) is a continuous, non-decreasing, and concave function of the traffic 
volume such that di{0) = 0; Ci = Ci{t) > 0 is the possibly time-varying maxi¬ 
mum flow capacity; = ai{t) G [0,1] is a possibly time-varying demand control 
parameter actuated via speed limit control on the non-source cells i G £ \ 7Z 
and ramp-metering on the sources i G (See Figure [^) On the other 

hand, the supply functions Si{xi) of every non-source cell i G £ \ TZ are as¬ 
sumed to be continuous, non-increasing, concave, and such that Si(0) > 0, with 
= infjxi > 0 : Si{xi) = 0} denoting cell i’s jam traffic volume. Conven¬ 
tionally, Si{xi) = - 1-00 at all sources i gTZ. 

®In the context of freeway networks, @-0 can be realized through appropriate setting 
of speed limits and ramp metering. In particular, for linear uncontrolled demand functions 
di{xi) = ViXi, formula § is equivalent to the modulation of the free-flow speed Vi = Viai, 
where Vi could be interpreted as the maximum possible speed due to, e.g., safety considerations 
(Cf., e.g., [H]). On the other hand, 0 corresponds to metering the maximum outflow from 
the onramp, which is its demand, by imposing a maximum value Ci = O-iCi (Cf., e.g., [22]). 
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Finally, we consider a nonnegative, possibly time-varying, £ x £ routing 
matrix R = R{t) satisfying the network topology constraints 


Rij =0, {i,j)e{£x£)\A, 


(5) 


and such that 

^ = 1 , tG£\S. ( 6 ) 

jes 

The matrix R is to be interpreted as describing the drivers’ route choices, with its 
entries Rij, sometimes referred to as turning ratios, representing the fractions of 
flow leaving cell i that is directed towards cell j. Equation ®> then guarantees 
that all the outflow from the non-sink cells is split among other cells in the 
network, while equation (i guarantees that the outflow from cell i is split among 
adjacent downstream cells only. Depending on the specific problem considered 
in this paper, the matrix R is to be considered as a control variable, as is the case 
for the Dynamic Traffic Assignment (DTA) problem, or an exogenous input, as 
is the case for the Freeway Network Control (FNC) problem. 

The traffic dynamics is then described by coupling the mass conservation 
laws 0 and 0 with the following functional dependence of the exogenous 
outflows and the cell-to-cell flows on the traffic volume and control parameters: 


— dk (,3^k , ^k^ , k ^ S , 

fij — RijZi , Zi — 'Yi (ii{xi,(y.i^ , (A j) ^ ; 

where for alH S f 


(7) 

( 8 ) 


'yf = sup <76 [ 0 , 1 ] : 7 - ™ax E Rhkdh{xh,oth) < Sk{xk) > (9) 

[_ {i',k)eA J 

Equations 0-0 generalize Daganzo’s cell-transmission model [S] by extending 
it to the case where junctions may have multiple incoming and outgoing cell^ 
and the cells’ demand and supply functions are allowed to be concave rather 
than piecewise linear functionslj Since fij = RijZi in every circumstances, equa¬ 
tions 0-0 enforce the FIFO^irst-In First-Out) constraint and are therefore 
amenable to the modeling of multi-source multi-sink transportation networks. 
For a given assignment 

Xi(0) = Xj , i € £, (10) 

of initial traffic volumes in the cells and dynamics of the inflows A(t), routing 
matrix R(t), and demand control parameters a(t), the evolution of the traffic 

^Daganzo [9] considers only the cases of nodes with single incoming link (diverge junction) 
or single outgoing link (merge junction). It is easily seen that, for diverge junctions §-§ 
prescribe that the outflow from a diverge junction always splits into the downstream cells 
according to the turning ratios, and the supply of a congested merge junction is allocated 
according to a proportional rule. This is a slight variation with respect to the original CTM 
model [^, in which merge is solved via a priority rule. 

^Daganzo’s original cell-transmission model [8] assumes a trapezoidal Fundamental Dia¬ 
gram, hence piecewise linear demand and affine supply functions. 
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volume vector x{t) for t > 0 is uniquely determined by equations Q, ([^, and 
0 -@. In particular, for every choice of and of X{t), R{t), and a{t) for 
t > 0 , the traffic volumes satisfy Xi{t) > 0 for every cell i G £ (since, when 
Xi = 0 , then Zi < di(xi, ai) < di{Q) = 0 ) and Xi{t) < 2 ;^“™ for every non-source 
cell i G £\TZ (since, when Xi = then < Si(xj'^"^) = 0). 

Within this setting, the DTA and the FNC problems can be cast as open- 
loop optimal control problems consisting in the minimization of the integral of 
a running cost ip{x, z) that is a function of vector of traffic volume and outflows 
over the time interval [0,T], where T > 0 is a given time horizon. We focus on 
running costs ipixjz) that are convex in (a:, z), nondecreasing in each entry x^, 
nonincreasing in each entry Zi, and such that 'ip{0,0) = 0. A relevant special 
case is when the cost function is separable, i.e., when 

'4’{x,z) ='^1pi{Xi,Zi) , ( 11 ) 


for convex costs 4’i{.Xi,Zi) of the traffic volume and outflow on the single cells 
i G £, with ^/>i(0, 0) = 0. This class of cost functions includes several standard 
choices [221 [I2] such as 

• Total Travel Time, for which ijji{xi,Zi) = xt. In fact, Xi{t)dt can be 
interpreted as the total time spent on cell i by all travelling vehicles in the 
interval [0,T]. 

• Total Travel Distance, for which 'il)i(xi,Zi) = —liZi, where ii stands for 
the length of cell i. In fact, £iZi{t)dt is the the distance travelled on 
cell i by all vehicles that exit from it in the interval [0, T], The minus sign 
implies that such a cost function is to be maximized, instead of minimized; 

• Total Delay, for which ipi{xi, Zi) = Xi — Zijvi, where Vi is the freeflow speed 
on cell i. In fact, the integral (xi(t) — Zi(t)/vi)dt represents the total 
additional time spent by vehicles on cell i with respect to the freeflow case 
Zi = ViXi during the interval [0,T]. 


Furthermore, any linear combination of such costs remains convex and might 
be used to trade-off between them. 

Given a single-sink transportation network, a running cost tp{x,z) as above, 
a finite time-horizon T > 0, initial cell traffic volumes x° > 0, and inflows 
Xi(t) > 0 for 0 < t < T at the source cells i G TZ, the DTA problem for the 
continuous-time cell-based dynamic traffic model of Section can be formulated 
as follows: 

rT 


min 

(y.{t) ,R{t): 


tf}{x{t), z{t))dt. 


( 12 ) 


Similarly, the FNC problem in a possibly multi-sink network can be formulated 
as follows, for a given running cost z), finite time-horizon T > 0 , initial cell 


traffic volumes > 0, exogenous inflows Xi{t) > 0 at the sources i S TZ and a 
routing matrix R{t) satisfying ([^-([^ for 0 < t < T: 


mm 




(13) 


Observe that the optimal control formulations (12) for the DTA and (131 


of the FNC problem differ because of the role of the routing matrix that is 
an endogenous optimization variable in the former, while it is exogenously de¬ 
termined in the latter, 
problems 


Both (12) and (131 are non-convex optimal control 

The 


which are hard to be either analyzed or numerically solved, 
source of non-convexity are the equations 0-® which are nonlinear equality 
constraints. In the next section, tight convex relaxation of the optimal control 


problems (12) and (131 will be discussed. 


3 Tight Convexifications of the DTA and the 
FNC problems 

In this section, we show how, in spite of their non-convexity, the DTA and 
the FNC problems (12) and (13) introduced in Section are amenable to be 
suitably reparametrized and relaxed in such a way that the obtained relaxations 
are both convex and tight, i.e., their optimal solution can be mapped back into 
a solution of the original optimal control problems with a proper choice of the 
demand control parameters (i.e., ramp metering and variable speed limits) and, 
for the DTA problem only, of the routing matrix. As will be clarified later, the 
approach is a generalization of the ones proposed in m and [T^ . 

We start by presenting below two optimal control problems that are convex, 
and are thus amenable to be analyzed and solved in a computationally efficient 
way. We then prove how an optimal solution of such convex problems can be 


mapped into one of the original problems (12) and (13). The convex optimal 


control problems that we introduce do not directly involve the control variables 
a{t) and R{t) as and ( [I^ , but rather the flow variables f{t), y{t), and z{t) 
along with the traffic volume vector x{t). Besides the law of conservation of 
mass ([^ and ([^ and the initial traffic volume equation (10), such variables are 
required to satisfy the nonnegativity constraints 


fij{t)>0, (j,j)eA, yi{t)>0, iGS, /Xi(t)=0 

and the supply and demand constraints 

yi{t) < Si{x^{t)) , z^{t) < mm{dt{xi{t)), CJ , 


i £ £ \ S, 
(14) 


i £ £. 


(15) 


The inequalities in ( |14[ ) enforce non-negativity of the cell-to-cell flows f tj and of 
the external outflows while the equality in (14) guarantees that outflows yi 


^Observe that, together with 0 and non-negativity of the exogenous inflows Aj, equation 
10 implies non-negativity of the cells’ inflows yi and outflows Zi as well. 
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towards the external world are possible only from the sink nodesj^ On the other 
hand, the inequalities in (15) capture the physical constraints on the cells: they 


guarantee that the total inflow yi in a cell does not exceed the supply Si{xi), 
and the total outflow Zi from a cell does not exceed the demand min{(ij(a;i), Ci}. 
Because of the assumption di(0) = 0 and non-negativity of the cell inflow y^, 
equation Q implies that the traffic volume Xi{t) remains nonnegative in time 
on every cell i G £. Analogously, for non-sink cells i G £ \ S, non-negativity of 
the outflow Zi and the fact that Si{xl°‘"^) = 0 imply that the traffic volume Xi{t) 
never exceeds the jam volume 

The first convex optimal control problem that we consider is a relaxation of 


the DTA (12) that reads 


mm / 

x. y,z , f.ij: Jq 

, |iot , , 


ip{x(t), z{t))dt. 


(16) 


By including an additional constraint ensuring that the outflow from cell i be 
split exactly as prescribed by an exogenous routing matrix R{t), i.e.. 


= R^j{t)z^{t), 


i G £, 


(17) 


we obtain a second convex optimal control problem 


min / 'tlj{x{t), z{t))dt. 


(18) 


that is a relaxation of the FNC problem (13). 


Convexity of the optimal control problems (fl^ and (fl^ means that, if 
(a;(°) (t),y(°^ it),z^°'> (t), (t),/r(°) (t)) and (a;(^) {t),z^^t), /(i) (t)) 

both satisfy the constraints in ( [I^ (respectively, in ([I^), then, for every /3 in 
[0,1], also {x^^\y''^\ z''^\ does, where x^p> = (1 — f3)x^°'> + /3x''^\ 

yW = (1 — j3)y^^'^ + Py^^^ and so on, and 

f ilj{x^^\t))dt < {1 — P) ( tlj{x^^\t))dt + P ( 'ijj{x^^\t))dt. (19) 

Jo Jo Jo 

In fact, it can be easily verified that concavity of the supply and demand func¬ 
tions implies that is a convex inequality constraint, while the remaining 
constraints (fTo|,([T4|, and ([Tt]) are all linear equalities, and convexity of 


the cost function ij}{x,z) implies (19). 


^Notice that, because of the assumption that Ai 7^ 0 only on source cells, and of the first line 
of ( |14[ l, it turns out that at most one between the two inflow terms Aj and fji appearing 
in the righthand side of the first equation in § can be positive: the exogenous inflow Aj for 
onramps i £ and the aggregate inflow from other cells fji for every i £ S\ 7 Z. Similarly, 
HI implies that only one of the outflow terms /ij and fij appearing in the righthand side 
of the second equation in § can be positive: the external outflow for off ramps i £ and 
the aggregate outflow towards other cells fij for all f G \ * 5 . 
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On the other hand, the fact that the optimal control problems (16 1 and (18 1 
are relaxations of and ( [T^ , respectively, is easily verifiable. Indeed, for 
every choice of demand control parameters a{t) in [ 0 , 1 ]^ and routing matrix 
R{t) satisfying (§ and (pL the traffic dynamics generated by Q-lH) and 0- 
( 10 ) necessarily satisfies (|14|), (15), and ( [T7| ), hence it is a feasible solution of 
the convex optimal control problems and That these relaxations are 
tight is implied by the following result, showing that for every feasible —hence, 
in particular, for the optimal— solution of (16) (of @) there exists a choice 
of demand control parameters a{t) and routing matrix R{t) (a choice of the 
demand control parameters a(t) for every exogenous routing matrix R(t)) such 
that equations 0-0 are satished, so that the solution is also feasible for ( [T^ 
(respectively, (13)). 


Proposition 1 Let Q = (V,£) be a network topology, a vector of initial 
traffic volumes, and Xi{t), for t > 0, exogenous inflows at the sources i € TZ. 
Then, 


(i) For any feasible solution {x{t),y{t), z{t), fi{t), f(t)) of the convex optimal 
control problem ( |16[ ), set the demand controls a{t) and controlled routing 
matrix R{t), for t G [0, T], as follows 


aflt) 


zflt) / dflxflt)) iGSXTl 

Zi{t)/C^, iGlZ. 


( 20 ) 


Rij (t) 


f,j{t)/zi{t) {i,j)GA 

0 {i,j)GSxS\A, 


( 21 ) 


with the convention that ai(t) = 1 if zflt) = dflxift)) = 0 on a non-source 
cell i G £\TZ, and that, if Zift) = 0, then Rijft) = \{k G £ : {i,k) G 
for all {i,j) G A. Then, R{t) satisfies the constraints 0-0 and x{t) 
satisfies the controlled dynamics 0 - 0 , so that {a{t),R{t)) is a feasible 
solution of the DTA problem 

Moreover, let R{t), t G [0,T], be an exogenous routing matrix satisfying 0 and 
0 . Then: 


(ii) For any feasible solution (x{f),y{f), z(t), p,lt), fff)) of the convex optimal 
control problem (18), set the demand controls aft), for t G [0,r], as in 
(20). Then, x(t) satisfies the controlled dynamics 0 - 0 , so that a{t) is 
a feasible solution of the FNC problem (13). 


Furthermore, in both cases (i)-(ii) above, the implemented trajectory is always 
in free-flow, i.e., zflf) = difxift),aiff)) for all t G [0,T]. 


The proof of Proposition [T] is provided in Appendix |A.1| See Figure for a 
graphical interpretation of the chosen demand controls afli) on the source and 
on the non-course cells, respectively. 
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Proposition provides a methodology to take any feasible solution of the 
convex optimal control problem (respectively, @) , and make it feasible 
with respect to the DTA problem (12) (respectively, to the FNC problem 
The convex optimal control problem ( |16| ) is a continuous-time version of the 
DTA formulation considered by m in discrete time and in the special case 
when the cost 

= ( 22 ) 

i^£ 

is the total traffic volume in the network and the demand functions di{xi) and 
supply functions Si(xi) are piecewise affine. In this case —and more in general, 
for linear, not necessarily identical, cost functions— (161 is a linear program. 


In our more general formulation, where the cost ip(x) is allowed to be a convex 
function of the traffic volume vector x and the demand and supply functions are 
concave, (16) is a convex (infinite-dimensional) program. The special case (|11[) 


where the cost is separable and convex is suitable for efficient solutions, e.g., 
based on distributed iterative algorithms [1]. On the other hand. Proposition]^ 
generalizes existing results, e.g., see [52], which are applicable only in specific 
scenarios, to convex costs, concave demand and supply functions, and arbitary 
network topologies. 


4 Analysis of the FNC problem 

In this section, we apply optimal control techniques to the study of the convex 


formulation of the FNC problem (18). Throughout, we focus on the case when 


the running cost il’ix) is a function of the traffic volume vector only. As a first 

step, it proves convenient to use the outflows Zi as the only control variables and 

substitute = RijZi for (i, j) G A, fik = Zk for k G S, and yi = Xi + 

for i G £. Then, the convex optimal control problem (18) can be reformulated 

as 

min / tp{x{t))dt 

Jo 

Xi{0) = x° , Xi{t) = X,{t)+ '^Rji{t)zj{t) - z^{t), iGS, 

j&£ 

0 < z^{t) < inm{d^{xi{t)),C^} , Rji{t)zj{t) < Si{xi{t)), iGS. 

3^£ 

(23) 


The Hamiltonian associated to the optimal control problem (231 is given by 

H{x, z, C) = i^ix) + '^q(^X^+Y^ RjiZj - z^ , (24) 


i£S 




where x is the state vector, z is the control vector, and ^ is the adjoint state 
vector. Upon introducing the notation 


— ^i{t) 'y RijCj{t) J 
je£ 


i G £, 
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the Hamiltonian can be rewritten as 


H{x, z, C) = tp{x) + QXi - Y 


KiZi 


(25) 


i^£ iGS 

Then, the Pontryagin maximum principle implies the following necessary con¬ 
dition: if {x*(t), z*(tj) is an optimal solution of (23l, then, for every t G [0,r], 


-Zi ■ 


z*{t)G argmin H{x*{t), zXit)) = argmax 

{Zi)iG£ iG£ 

0 < Zi < Ci 0 < Zi < Ci 

Zi<di{x*{t)) 0 < Zi < d^{x*) 

y ' ^jjZj ^ £iiXi{t)) y ( ddjjZj < Si(x^) 

j&£ je£ 

(26) 

For a given value of the optimal traffic volume vector x*{t) and of the adjoint 
state vector ^(t) (hence of K{t)), the optimization in rightmost side of ( [^ is a 
linear problem in the variables Zi. If one denotes by Vi, and Xi the multipliers 
associated to, respectively, the demand, the supply, and the capacity constraints, 
then the dual problem to one in the rightmost side of (261 can be written as 


G argmin +’^^Si{x*{t)) + XiCt) , 


e, > 0 , Z/, > 0 , > 0 

+ Xi + X! ^ 

3&£ 

and the adjoint dynamical equations read 


(27) 


at) = -^ax*m ++ v*{t)s'M{t)) , 


iG8, (28) 


where dt{xi) and st(xi) are the derivatives of the demand and, respectively, of 
the supply functions, with transversality condition 


0 (r) = o, 


i G £. 


(29) 


If every node v in the network is either an ordinary junction (single incoming and 
outgoing cell), a merge junction (multiple incoming and single outgoing cells), 
or a diverge junction (single incoming and multiple outgoing cells), then it is 
convenient to regroup the addends of the summation appearing in the rightmost 
side of (261 so that the necessary condition for optimality can be separated into 
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decoupled local linear programs 


Ki{t)zi , Vn G V. 


G argmax 

0<Zi<Ci 

Zi < d^ix*it)) 

RijZi < Sj{x*{t)), \/j■■ <jj=v 

i: Ti—v 

(30) 

In the special case where the cost is the total traffic volume, the increasing 
parts of the demand functions are linear and the supply functions are affine, 
all with identical slopes, equations (28)—(30) imply the following result on the 
structure of the optimal control. 


Proposition 2 Let the cost function be as in (22), the increasing parts of the 
demand functions be di{xi) = uiXi, and the supply functions have the form 
Si{xi) = 6i — LoXi for some 0^ > 0 and w > 0. Let R be a routing matrix 
that is constant in time, and let {x* ft), z*(t)) be an optimal solution of the 
corresponding FNC. Then, for every non-sink cell i G S \ TZo, 

(i) if V = Ti is an ordinary junction with downstream cell j (with aj = v), 
then 

z*{t) =miii{d^{x*{t)),Sj{x*ft))}-, (31) 

(ii) if V = Ti is a diverging junction, then 

zfit) = if ft)d^{x*{t)) 


if W = sup {7 G [0,1] : jR^kdi{x*{t)) < Sk{xl{t)), Vk G S} ; 


(32) 


(Hi) if V = Ti is a merging junction with downstream cell j (with aj = Ti), then 

{zk{t))h-.Th=v G argmax ^ Kh{t)zh{t). (33) 

0 < Zk < dkixlit)) 

Sj{x*{t)) 

hes- 

Th—V 


Proposition [^implies that at diverging junctions of a transportation network 
with total traffic volume as cost, affine supply functions and linear demand with 


identical slope, the optimal solution satisfies (321 which coincides with the FIFO 


diverge rule of Daganzo’s cell transmission model. Similarly, for merging junc¬ 


tions with two upstream cells h and i and downstream cell j, (33) is equivalent 
to Daganzo’s priority rule 


/*■ = mid{d*(a:-),Si(a;*) - dh{xl),p^Sj{x*)} 
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and 


fhj = inid{dh{xl),Sj{x*) - d^{x*),phSj{x*j)} , 

where mid{a, 6 , c} denotes the median and the priority parameters ph and pi 
are such that: Pi = 1 and ph = 0 if Ki > Pi = 0 and ph — f if Ki < Kh', and 
Pi = 1 — Ph is arbitrary in [0,1] if Ki = k/i 

Proposition and the discussion following it imply that, if the demand and 
supply functions are linear with identical slope, the optimal solution of the FNC 
problem with total traffic volume as a cost, is achieved by the CTM with 
FIFO rule at the diverge junctions with no additional control required, i.e., 
ai{t) = 1, at all the cells i immediately upstream of a diverge or an ordinary 
junction, while additional control (i.e., ai{t) < 1 ) might be useful on the cells 
immediately upstream the merge junctions. Observe that, for costs different 
than the total traffic volume, additional control might be beneficial at any cell, 
being it immediately upstream a merge, diverge, or ordinary junction. 


5 Robustness analysis 


In Section]^ we discussed how the DTA (12 1 and the FNC (13) can be cast 
as convex optimization problems (16) and (18) whose solutions can be mapped 
back into optimal demand control parameters a(t) and routing matrix R(t) as 
explained in Proposition]^ In order to solve (16) and (18) and then find com¬ 
pute the optimal a{t) and i?(t), one needs precise information about the input 
parameters, i.e., the initial traffic volume x^ and exogenous inflow vector A(t) 
over the planning horizon [OiT]. However, in practice, information about these 
quantities inevitably involves uncertainties. Therefore, a reasonable strategy is 
to (i) compute the optimal solution of or ( |13[ ) for nominal values of the 
exogenous inflows and initial traffic volume, and (ii) compute the corresponding 
nominal control inputs using Proposition]^ and then expect the trajectory un¬ 
der actual parameter values and the nominal control inputs to be close enough 
to the nominal trajectory. In this section, we provide formal guarantees on the 
robustness of a general system trajectory to perturbations in the initial traffic 
volume and the exogenous inflows, while maintaining the same control input]^ 
When specialized to the system trajectory corresponding to the optimal solution 
of the DTA (12) or the FNC (13) under nominal values of the initial traffic vol¬ 
ume and exogenous inflows, this gives the desired result on robustness analysis 
with respect to uncertainties in the input parameters. 

We shall use the notational convention that x^ and A(<) denote the nominal 
values of the initial traffic volume and the exogenous inflows, while x^ and 
X{t) denote the perturbed values of these parameters. Similarly x{t) and x{t) 
will, respectively, denote the nominal and perturbed trajectories, under the 
same open-loop control inputs a{t) and R(t). Our robustness analysis will 


®In fact, when Ki = k^, the proportional rule fij = di{x*) min{l, Sj{x*)/{di{x*) + 
dh{x’^))}, fhj = dh{xDm:m{l,Sj{x*-)/{di{x*) + dh{xD)} also satisfies l |^ . 

®We recall that the control inputs are open loop, and not in feedback form. 
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provide bounds on perturbations in the state trajectory due to perturbations in 
the inflow A — A and the perturbations in initial traffic volume — x^. The 
bounds derived in the general setting, when specialized to the case where x{t) 
corresponds to the optimal solution of the DTA (12) or the FNC (13), and where 
a{t), and R{t) are derived from x{t) according to Propositionwill give the 
desired robustness analysis for optimal solutions. 

Our technique relies on leveraging a certain monotonicity property of the dy¬ 
namical system underlying the dynamics in -(|l0|). Monotone systems 

are dynamical systems whose trajectories preserve the partial order between 
initial traffic volumes and external inputs. Specifically, for monotone systems, 
trajectories with initial traffic volumea^a;° < and inflows X{t) < \{t) for all 
t e [0,T], are such that x{t) < x{t) for all t e [0,T]. For given a{t) and R{t), 
let us rewrite the dynamics in compactly 


as: 


ii = Xi(t) + g^{x, a, R) 


(34) 


A standard result in dynamical systems theory, known as Kamke’s theorem 


usi Theorem 1.2], then implies that ( |34[ ) is monotone in a certain domain R C 
if and only if 

fin ■ 

(35) 


^{x,a,R) > 0 , 

dXj 


for almost every x G V, at all t G [0,T]. It has been recognized, e.g., see [S], 
that, for a given a{t) and R{t), the dynamical system (341 is monotone, if x{t) 
in the free-flow region, i.e., it satisfies x G where Rijdi{xi, at) < Sj(xj) 
for all j G £ for all t G [0,r]. However, the monotone property is not satisfied 
in the region where there exists at least one cell k outgoing from a diverge 
junction (i.e., a cell j such that there are other cells k ^ j with ak = crj) that 
is congested, i.e., Rijdi(xi,ai) > Sj{xj). In fact, in such a configuration, 

monotonicity is lost since an increase of the traffic volume Xj in a congested cell 
j has the effect of reducing its supply Sj{xj), and thus the coefficient (as 
defined in 0) of any upstream cell i such that {i,j) G A. In turn, this implies 
a reduction of the flow fik from any such cell i to any other downstream cell k 
such that {i,k) S A, hence a reduction in the inflow on cell k. 

Given an n-dimensional vector x, recall that its £i norm is given by ||a:||i := 
Er=i l®*l- For example, ||i° - a;°||i = “ *Wlli = 

~ for all t G [0,rj. We let ||A — , 


Im :=/J|lA(s)-A(s)||ids. 


Proposition 3 If the perturbations 
small, then 


and IIA — A||i,t are sufficiently 


||i(t)-x(t)||i < ||i°-a:°||i + ||A-A|li,t, VtG[0,T]. (36) 

The sufficiently small condition in Proposition is to ensure that the per¬ 
turbed trajectories remain in free-flow, and hence the system is monotone along 
them. The bound provided by Proposition]^ has the desirable property that the 

^'^Here by a < 6 for vectors in we mean at < bi for all i £ £. 


16 







right hand side of (361 goes to zero as the magnitude of perturbations go to zero. 


However, the bound is conservative for non-negligible values of perturbations. 
The next result provides a tighter bound in this latter case. In preparation for 
it, let ||Ai - AjIIoo := sup^gjo-^] |Aj(t) - Aj(t)|, i G Tl, and 

Ai := sup {Ai(t) + ||Ai-Ai||oo}, Aj := max{0, inf {Ai(t)-||A — AiHoo}}, *G 

tG[0.T] ‘G[0,T] 


:=x° + |x°-x°| 


x° := max{0,x° - |x° - x°|}, i&S 


Note that the extreme values of perturbed inflows, A and A, are constant. 
The next bound is expressed in terms of equilibrium traffic volumes at these 
extreme values, denoted as x®‘5(A) and x®'i(A). Existence and uniqueness of 
these equilibria are discussed in [5]. In particular, if an equilibrium exists, it is 
unique. 

Proposition 4 —x°||i and ||A —A||i^t sujficiently small, and 

and x®®(A) exist, then 


\x{t) - x{t)\\i < ||x®«(A) - x®«(A)||] 


||x°-x°||i- 


min {||x®«(A) - eili + ||a:^«(A) - e||i} Vt G [0,r]. 

(37) 

Again, in Proposition the smallness of the perturbations is required to 
ensure that the perturbed trajectories remain in the free-flow region and that 
the corresponding extreme values of inflows admit equilibria. 


Remark 1 The optimal solution x(t) to the DTA (12) or the FNC (13) under 
the control policies (a, R) designed in Proposition is in the free flow region. 
Hence, Propositions and are applicable when the nominal trajectory x{t) 
corresponds to the evolution of the system trajectory under the control policies 
designed in Proposition^ 


Note that the right hand side of (37) does not go to zero as perturbations 


go to zero. However, for non-negligible perturbations, it gives less conservative 


upper bound than (|36|). Therefore, we propose to use the bound equal to the 

and 


minimum of the right hand sides of ( |36[ ) and ( |37[ ). Indeed, this is the upper 
bound we use for comparison with simulation results. 

The perturbations for which Propositions|^and|^are not applicable typically 
result in an overload condition, where the perturbed exogenous inflow exceeds 
network capacity. In such a scenario, the traffic volume on the sources grow un¬ 
bounded. In our previous work [3l Proposition 2], we characterized this growth 
rate for a related dynamical network flow model. Inspired by our previous work, 
we suggest a simple upper bound in the overload regime, for networks with a 
single source, e.g., as shown in Figure]^ and with constant values for nominal 
and perturbed exogenous inflow. The suggested upper bound will be compared 
in the next Section against simulation results. Let A be the supremum of all 
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perturbed exogenous inflows under which the perturbed trajectories remain in 
free-flow, and hence for which Propositions!^ andare applicable. Let x denote 
the perturbed trajectory under perturbed inflow A. Triangle inequality implies 
that 

||i(t) - a;(t)||i < p(t) - a;(t)||i + ||S(t) - i(t)||i Vt>0. (38) 

The first term in (38) can be upper bounded using Propositions and In 
the second term, the total traffic volume on sources grows unbounded at a rate 
which is expected to be equal to ||A — A||i = ~ -^*1: i-®-; 


t—¥00 ^ 


(39) 


Since (391 gives a reasonable upper bound for large t, considering the relatively 
large time horizon for the simulations, we shall use 

||i(f)-i(t)||i < ||A-A||it (40) 

for all t. In summary, the upper bound for the overload regime is obtained from 
([38|) and (|40|. 


The upper bounds in (361, (371 and (401 are to be contrasted with the bounds 
obtained from standard sensitivity analysis of ordinary differential equations: 


^ e-^»*||A(s) - A(s)||i ds + ||i(0) - x(0)||ij , t G [0,r] 

(41) 

for any Lg > 0 such that, for all t G [0, T], 
llg(x,a(t), E(t))-g(x, a(t), < Lg\\x-x\\i , Vx,a; S n,e£[0,xf“] (42) 


The derivation of (41) is provided in Section A.5 A particular choice of Lg 
is also provided in Lemmaj^in Section A.5 If ||A(t) — A(t)||i is constant over 


[0,T], then (41) becomes 

pLgt - 

\~x{t)-m\\i< 


||A-A|li+e^«‘||i(0)-x(0)|li, <e[0,T] (43) 


The bounds that we have obtained are very easy to compute and yet pro¬ 
vide tight upper limits on the cost displacement as a function of the inflow 
uncertainty, as we illustrate via numerical examples in Section [6.2| 


5.1 Extensions to Non-FIFO Models 

The proportional rule in (§ for traffic flow under congestion at nodes with mul¬ 
tiple outgoing links is often also referred to as the FIFO rule. While such FIFO 
rules are a natural framework for multi-origin multi-destination transportation 
networks, it has been pointed out by several authors, e.g., see |T0], that it cor¬ 
responds to a rather conservative behavioral model, in which not only drivers 
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never change their routing choice, thus blindly queuing up even in presence of al¬ 
ternative routes to the same destination, but more importantly it does not take 
into account presence of multiple lanes for multiple maneuvres at junctions. For 
example, a congested offramp on a freeway would slow down and possibly block 
the flow of vehicles on the main line, which is not always realistic. Recently, 
e.g., see ini, there has also been interest in non-FIFO rules, for which § is 
replaced by 




7 = sup < 7 € 


[oa] : 7 • X! ^hjdhixh) < Sjixj) 


(44) 


hee 


The main difference of the non-FIFO model with respect to the FIFO one is 
that congestion in one of the outgoing cells does not influence the flow towards 
other outgoing cells. In the previous example, while the congested offramp 
forces vehicles that would like to take it to stop in the freeway, those that want 
to continue on the main line are free to do so, if the downstream supply of the 
main line is sufficient. The non-FIFO model is best suited for single-origin single¬ 
destination network (in which the actual path of a vehicle does not matter), 
but could possibly exhibit unrealistic behavior in multi-origin multi-destination 
networks. A possibly more realistic model may involve a combination of FIFO 
and non-FIFO, in which a fraction of drivers can change path, while another 
fraction cannot - for example, private cars can deviate from their path to chose 
a more convenient one, while buses or trucks have prescribed paths to follow. 

Traffic systems under non-FIFO models can be shown to be monotone almost 
everywhere in 27 = e.g., see m- This is because an increase in traffic 

volume in an outgoing link does not decrease the flow towards other outgoing 
links at the same junction, as thoroughly discussed in [2] . As such, their stability 
properties and their set of equilibria are substantially different from a traffic 
system under FIFO model m- However, several results that we stated for 
FIFO models can be naturally extended to non-FIFO models. In fact, the 
control design derived in Proposition holds true for traffic dynamics Q-©, 
0 , (§ and ( |44[ ) with no difference, since the resulting optimal trajectory is in 
free-flow, where = 7 ^ and where thus FIFO and non-FIFO models coincide. 
Moreover, the robustness analysis of Section can be extended, for a given 
a(t) and R{t), as follows. Since the system is monotone everywhere in 27 = 


the robustness bound in (36) holds true independent of the magnitude of 

In our previous work [ 


perturbations ||x° — x°||i and ||A — A|| 


x“ — a:“Hi ana ||A — in our previous worK we have 

shown that, under the non-FIFO model, in general, there exists a manifold of 
equilibria for constant exogenous inflows A. Nevertheless, the upper bound in 
( |37[ ) is valid for any choice of equilibrium. Numerical simulations reported in 
Figure and [^suggest that the non-FIFO dynamics tend to be more robust to 
perturbations than the FIFO one. 
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Figure 4: The network used in the numerical study. 


6 Numerical studies 

6.1 Simulations Comparing Various DTA Formulations 

In this section, we compare the optimal solutions obtained from the two pro¬ 
posed relaxations of the DTA problem and of the FNC problem ( [T^ , with 
an uncontrolled traffic simulation using FIFO model, following [^, but with 
proportional merge rule as in |22[ [6], as described in Section We remark 
once again that the optimal solutions are not meant to represent, on their own, 
traffic simulations, but rather optimal trajectories, which can be made feasible 
with respect to the traffic dynamics via the control signals that are designed in 
Proposition 

We first provide the details of the simulation setting that are common for all 
our numerical studies. We solve the convex optimization problems in MatLab 
using the Convex Programming package cvx mm- We use the single-source 
single-sink network described in [27] and shown in Figure]^ for our simulations. 
For implementation, we discretize the continuous formulation according to stan¬ 
dard practices in Cell Transmission Models. Time is slotted with sampling time 
r = 10 seconds. In all the cells, demand and supply functions are piecewise 
affine, with 




Si{xi, t) — min • 


/ lam \ 

[xi - Xi) 




where Vi, Wi, Ciit) (t = 1,... ,n), Li and are the free-flow speed, the wave 
speed, the capacity, the length and the jam traffic volume on cell i, respectively. 
Recall that the actual demand is the saturated mm{di{xi),Ci{t)}. The values 
of these parameters, along with number of lanes and length of cells, are specified 
in Table The units of all parameters, as well as of inflows, provided below, 
are chosen in such a way that physical consistency is ensured. In addition, 
with the chosen parameters, a vehicle travels along an entire cell in exactly one 
time slot at maximum speed, which in the considered scenario is the free-flow 
speed V. Therefore, the Courant-Friedrichs-Levy condition < 1, which 

is necessary for numerical stability, is satisfied. 

Vehicles enter the network from cell 1 at rate Xi{t). The setup in this Section 
and Section[6.3| differs from the setup in Section[6.2|in terms of the time-varying 
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Parameter 

Free-flow speed Vi, wave speed Wi 

Length of cell 

Capacity Ci 

Number of lanes ii 

Jam volume 


Value 

50 feet / sec 
500 feet 

Mi veh/r (except 4) 

2 for z = 1,2,9,10; 1 otherwise 
IMi veh 


Table 1: Cell parameters. 


values of Xi(t) and capacities on the links. In this Section and in Section 6.3 
we consider a setting in whicl|^Ai(l) = 8, Ai(2) = 16, Ai(3) = 8 and Xi{t) = 0 
for t > 4, with a time horizon of T = 25 steps, and in which the capacity in 
the cells is constant except on cell 4, where a bottleneck is simulated by setting 
€ 4 ( 1 ) = 6 veh/r for t ^ 5,6,7,8, (74(5) = € 4 ( 6 ) = 0 veh/r, (74(7) = ( 74 ( 8 ) = 3 
veh/r. Exogenous turning ratios are as follows: R 23 = 2/3, R 25 = 1/13 and 
i?34 = 2/3, i?36 = 1/3, the others being trivial. 

For these values of inflow and link capacities, we compute the optimal tra¬ 
jectories and costs for the DTA and FNC problems and the costs associated 
with a system evolving under FIFO traffic dynamics with proportional merge 
rule. For brevity, we refer to the last model as FIFO. The initial traffic vol¬ 
ume for each case was x(0) = 0. The results for the total traffic volume cost 
4'7 )(x) = reported in Table whereas the results for the 

quadratic cost, 4'^^)(x) = reported in Table The corre¬ 

sponding trajectories for volume of vehicles, Xi{t), for a few representative cells 
are shown in Figuresand Figuresfor linear and quadratic cost, respectively. 

As expected, the DTA, being the least constrained, gives a cost that is 
smaller than the FNC scheme, and both these optimal solutions yield a total 
cost that is no larger than the cost computed from the traffic simulation under 
FIFO traffic dynamics. This favorable comparison between the optimal costs 
for the DTA variants and the FIFO model served as a motivation to investigate 
the feasibility of optimal DTA solutions with respect to FIFO traffic dynamics, 
which we addressed in Proposition 

Interestingly, for the linear cost criterion, the FNC optimal cost coincides 
exactly with the FIFO case, confirming Proposition which shows that this 
no mere coincidence, and that one can identify a class of settings for which this 
property can be proven to be true. Finally, we refer to the cost comparison for 
the quadratic cost criterion in Table to emphasize that this property does not 
hold true in general. 


Inflows are normalized by the onramp length as the dynamics involve the traffic volume 
of vehicles. 
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Figure 5: Trajectories of the number of vehicles on cells 1, 2,3,4 for the system 
under FIFO traffic dynamics and for the optimal solutions corresponding to the 
two variants of DTA, for linear cost. 


Scheme 

Cost 

FIFO 

281.6 

DTA 

246 

FNC 

281.6 


Table 2: Comparison between optimal cost for the two DTA variants, and the 
cost for the system under FIFO traffic dynamics, for the linear cost criterion. 


Scheme 

Cost 

FIFO 

1930.5 

DTA 

1393.5 

FNC 

1595.7 


Table 3: Comparison between optimal cost for the two DTA variants, and the 
cost for the system under FIFO traffic dynamics, for the quadratic cost criterion. 
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Figure 6: Trajectories of the number of vehicles on cells 1, 2,3,4 for the system 
under FIFO traffic dynamics and for the optimal solutions corresponding to the 
two DTA variants, for quadratic cost. 


6.2 Robustness bounds 


In this second numerical study we compare the robustness bounds obtained in 
Propositions and 1^ and in (38) and (40), with simulation results. We run 
simulations for a constant inflow scenario where the nominal value of the inflow 
is Xi{t) = 5 for alH = 0,1,...,T, T = 200, and in which the link capacities are 
as in Section 6.1 except that C 4 {t) = 6 veh/r for all t, i.e., we do not consider 
any bottleneck in this section. 

For these parameters, and for the total traffic volume cost 
we solve the FNC with initial traffic volume a;(0) = 0, and compute the cor¬ 
responding controls, as given by Proposition We then compute the cost for 
the FIFO and non-FIFO dynamic traffic models under these controls, but under 
perturbed inflow Xi{t) = Xi(t) + AA and, for simplicity, zero uncertainty in the 
initial traffic volume (namely, x(0) = i(0)). We emphasize that the controls do 
not change with AA. Let the optimal trajectory, i.e., the solution of FNC, be 
denoted as x*{t), and the trajectory under perturbed inflow A A be denoted as 
(Note that x°{t) = x*{t)). We compute the resulting perturbation in 
cost A'I'(AA) = i^) ~ various values of AA in [0,3]. 

We compared our bounds both for a system under FIFO traffic dynamics 
and for a system under non-FIFO traffic dynamics, and the results are shown 
in Figure and in Figure respectively. 

The solid line plots for both FIFO and non-FIFO traffic dynamics in are 
piece-wise affine, with the transition point corresponding to network capacity. 
This transition point is at AA « 0.8 and AA « 2.8 for the FIFO and non-FIFO 
dynamics, respectively. In fact, when AA is less than the transition point, the 
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Figure 7: Comparison of the perturbation in cost due to perturbation in inflow 
under FIFO traffic dynamics obtained from simulations (solid line), correspond¬ 
ing bounds from Propositions!^ and(38) and ( [40| (dashed lines), and (on the 
right panel) from Sensitivity analysis (43). 




0 0.5 


1.5 2 2.5 3 0 0.5 

AA 


1.5 2 2.5 3 

AX 


Figure 8: Comparison of the perturbation in cost due to perturbation in inflow 
under non-FIFO traffic dynamics obtained from simulations (solid line), corre¬ 
sponding bounds from Propositionsand(38) and ( [40| (dashed lines), and 
(on the right panel) from Sensitivity analysis (43). 
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trajectories (t) are found to reach equilibrium, thereby allowing application 
of Proposition!^ Moreover, in this regime, the trajectories are in the free-flow 
region. When AA exceeds the transition point, then the trajectories are in 
overload, and the traffic volumes on onramps grow unbounded, which results 
in a higher slope in the perturbation plot in both Figure and Figure For 
AA < 0.8 (corresponding to the transition point for the FIFO traffic dynamcis), 
the perturbed trajectories under both dynamics are identical. That is, the 
better robustness property of non-FIFO models can be attributed to a higher 
network capacity, which in turn can be attributed to flexibility in routing of 
traffic flow in congestion. These features are well captured by the analytical 
bounds computed by Proposition and (when an equilibrium exists) and ( [M| 
and (401 (in overload), which are shown as dashed lines in Figure]^ and Figure 
In general, the perturbed trajectories under controls derived from DTA or FNC 
optimal solutions show qualitatively similar behavior, with the difference being 
in the locations of transition points. 

Finally, the bound resulting from Sensitivity anaylsis and (431 is shown only 
on the right panel and in log-scale. Differently from the robustness bounds that 
we derive, the order of magnitude of the error provided by Sensitivity analysis 
is not correct, providing over-conservative upper bounds. The reason is that 
while our robustness bounds depend on the time horizon T in a linear way, the 
Sensitivity analysis provides a worst-case bound which grows exponentially with 
T. 


6.3 Robustness-Performance Tradeoff nnder Reduced Fea¬ 
sible Set for DTA 

From the analysis proposed in Section is can be noticed that the network con¬ 
trolled via Proposition exhibits less robustness to inflow perturbations when 
the perturbed trajectories are in the congestion region. In this last numerical 
study we investigate the possibility of scaling down the supply functions in the 
DTA formulations, whose optimal solution is used to set the controls. Such a 
modification introduces more slack, thereby increasing the magnitude of inflow 
perturbations beyond which the trajectories enter the congestion region, and as 
such increases the network robustness. Naturally, such a modification comes at 
the expense of increase in the value of optimal cost under zero perturbation. 
The objective in this section is to study this tradeoff. 

Specifically, we replace the first inequality in ( |15[ ) with 

Viit) < s,{xi){l - e) ( fl^ ) 

where 0 < £ < 1 is a tunable parameter. In this section, we focus on the controls 
computed from optimal solutions to FNC; the results are qualitatively similar 
for DTA. 

We consider again the total traffic volume cost il^ix) = J2iee 
scenario described in Section |6.1| For a given e > 0, similar to Section 
we computed the control values from the optimal solution to FNC with initial 
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traffic volume x(0) = 0, and then simulated the system with several values 
of perturbations to inflow AA. The simulations were repeated for a range of 
values of the parameter e. Let (t) denote the perturbed trajectory. For 

each combination of e and AA, we computed the perturbed cost: 4'(e, AA) := 

relatively short time horizon setup of this 
scenario, the perturbed cost does not necessarily reflect the congestion effects 
in the perturbed system. Therefore, we additionally record congestion factor: 
7(e, AA) = minj_t where the minimum is taken over v &V and 

all times, and where 7_,- is defined in Clearly, 7(2, AA) < 1, where equality 
implies that the perturbed trajectory is in free-flow all the time, whereas strict 
inequality implies congestion. The smaller the value of 7(e, AA) is, the more 
congested some cells are. 

The results are shown in Fig. As expected, for every AA, the perturbed 
cost is an increasing function of e. Concerning congestion factor, 7(e,0) = 1 
for all e, as the optimal trajectory is always in free-flow, while for high AA, the 
lower £ is, the more the system is prone to congestion. At the two extremes, 
for e = 0 a small inflow increase yields congestion; for e = 0.5, the perturbed 
solutions are in free-flow even for relatively high values of AA. 



1 

AA 



Figure 9: Impact of reduced feasible set and inflow perturbation. Left panel: 
cost '!'(£, AA) as a function of A A for different e (colored lines). Right panel: 
congestion factor 7(e, AA) as a function of AA for different e. 


7 Conclusion 

We consider two variants of cell-based continuous time System Optimum Dy¬ 
namic Traffic Assignment, which we refer to as DTA problem, in which turning 
ratios are a possible control input, and FNC problem, in which turning ratios 
are exogenously assigned. Within these formulations we study relaxations of 
traffic dynamics, where total inflow into and total outflow from the cells are 
upper bounded by supply and demand respectively, and which are shown to 
yield convex problems. We also design open-loop variable speed limits, ramp 
metering and routing controllers that ensure feasibility of the optimal solutions. 
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This significantly expands known results in terms of relationship between com¬ 
putationally efficient DTA formulations and feasibility of their optimal solutions 
with respect to traffic dynamics. Furthermore, using monotonicity arguments, 
we derive bounds on the cost increase due to perturbations of initial conditions 
and external inflows under the designed open loop controllers. The proposed 
methodologies are illustrated via extensive simulation results. 

Future research direction includes leveraging necessary and sufficient condi¬ 
tions along the lines of Section and [lain] to possibly speed up numerical 
solutions to various variants of DTA. We also plan to develop distributed algo¬ 
rithms for solving the DTA and FNC problems, along the lines of our preliminary 
work in [T]. 
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A Proofs 


A.l Proof of Proposition 

(i) Let {x{t), y{t), z{t), fi{t), f{t)) be a feasible solution of the convex optimal 
control problem (TT^. At any time t > 0, the demand constraints in (151 


imply that Zi(t) < Ci. It follows that, with the choice of demand control 
parameters in ([^, for every non-source cell i G £\Tl, 


Zi{t) = min{Zi(t), Ci} = \m\i{ai(t)di{xi{t)),Ci} = di{xi(f),ai(t)). 


For sink cells k G S, the above readily implies that yk{t) = z^ft) = 
dk{xk{t),ak{t)), so that 0 holds true. On the other hand, for every 
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source cell i £ TZ, the choice of demand control parameters in (20) and 
the demand constraint Zi(t) < di(t) in (15) imply that 

Zi{t) = mm{di{x^{t)),Zt{t)} = min{(ii(a;i(t)), ai{t)Ci} = ~di{xi{t), ai{t)). 

Then, Zi{t) = di{xi{t),ai{t)) for every cell i £ E and the choice of the 
routing matrix © implies that 

/y (t) = Rij{t)zi{t) = Rij{t)di{xi{t),a^{t)) , 

for every pair of adjacent cells (i, j) G A. Therefore, for every cell j £ £, 

R^jit)di{x^{t),a^{t)) = /uW = ^ Sj{xj{t)), (45) 


where the inequality is implied by the supply constraint in (15). It follows 
from § that yf (t) = 1 for every cell i £ £ and this implies that 


yf {t)R,j(t)d,{xi{t),a,{t)) = , {i,j) £ A, 


(46) 


so that (|^ is satisfied. Finally, the choice of the controlled routing matrix 
(21) can be readily verified to satisfy (H and §. Hence, for every feasible 
solution {x{t),y{t), z{t), f{t)) of the convex optimal control problem 

(16), the choices (20) and (21) of the demand control parameters a{t) and 
of the controlled routing matrix R{t) satisfy the constraints , so 

that {a{t)jR{t)) is a feasible solution of the DTA problem Q- 

(ii) Let {x{t), y{t), z{t), f{t)) be a feasible solution of the convex optimal 

control problem ( |I^ . Then, for every time t > 0, the choice of the demand 
control parameters (20) along with the additional constraint © imply 
that 

/y (t) = Rij{t)zi{t) = Rij{t)di{x,{t),a,it)) , 
for every pair of adjacent cells (i, j) £ A. It follows that 


Rijdiixiit)^ (Xiity) — ^ ^ ; 


j & £ ■. 


i^S 


where the inequality is implied by the supply constraint in (15). From this, 
follows as in (i), so that 


is satisfied. On the other hand, one has 
that = Zk{t) = dk{xk{t),ak{t)) for every sink cell k £ S, so that also 
0 is satisfied. Hence, for every feasible solution {x{t), y{t), z{t), y,{t), f{t)) 
of the convex optimal control problem ( |l8| ) , the choice (20) of the demand 
control parameters a{t) satisfie the constraints 0-0. so that a{t) is a 
feasible solution of the FNC problem (13). 


The fact that x{t) is in free-flow under the designed controls follows from (45) 
in both cases (i) and (ii). □ 
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A.2 Proof of Proposition 

If II = (ji is a merging junction with downstream cell j, then we immediately 
get that 

G argmax 

o<z^< a 

Zi < d^{x*{t)) 

'^Zi< Sj{x*{t)) 
i£S: 

Ti—V 

On the other hand, if is either an ordinary or a diverging junction, then 

(47) 


a minimizer in (27) necessarily satisfies 


X* (t) = 0 , (t) + (t) = max{0, K,{t)} . 

Indeed, reducing Xi > 0, > 0, and Vj > 0 for all j such that aj = Ti until 

Xi = 0 and + ~ niaxjO, Ki} reduces the cost 


Y + XiCi) 

i£S 


without violating any of the constraints in (27). (The fact that is either 
an ordinary or a diverging junction implies that such variables Vj appear in 
no constraint Xh + ih + Sjefother than for h = i.) Then, by 
combining (28) and ([47| we get that 


(48) 


Ci ^ Cj 




= -1 + ^ i?,,- + cu ^ + er - ^ 


je£ 

< wmaxjO, Ki}, 


je£ 


je£ 


where the last step follows from the constraints = Ij t'i > 0, and 

— 0- It then follows from (481 and (29) that 

n^{t)>Q, tG[0,T]. (49) 


Indeed, if Ki{t) < 0 for some t € [0,T], then (48) implies that ki{t) = 0, hence 
Ki(t') = Ki{t) < 0 for all t' > t, which would contradict (29). It follows that 

argmax KiZi = max Zi. (50) 

0 < z* < Q 0<z^<Ci 

Zi<di{x*{t)) Zi<d,{x*{t)) 

RijZi < Sj{x*{t)) RijZi < Sj{x*{t)) 

If Ti is an ordinary junction, then ( [50| reduces to If Ti is a diverging 

junction, then it coincides with (32). 
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A.3 Proof of Proposition 

The proof is based on application of contraction principles, developed in our 
previous work [3], for system trajectories in the monotone region. First note 
that 


dt 


\\x{t) - x{t)\\i =^sgn(ij(t) - x,{t)) (A,(t) 


A.(<)) 


+ ^ sgn {xi{t) - Xi{t)) {gi{x, a, R) 

iGE 

<l|A(t)-A(t)||i 


+ X! - Xi{t)) {gi{x, a, R) 

iGS 


g^{x,a,R) 


g,{x,a,R)) 

(51) 


For sufficiently small perturbations, x{t) and x{t) remain in the free-flow 
region for all t G [0,r]. Therefore, the traffic flow dynamics is monotone, hence 
[3l Lemma 1] implies that the second term in the RHS of (51) is non-positive. 
Therefore, 


dt 


\x{t) - x{t)\\i < ||A(<) - A(t)||i Vt e [0,T] 


which, upon integration, gives (36). 


A.4 Proof of Proposition 

Let A) denote the state of the network at time t starting from initial traf¬ 

fic volume x'^ and under inflow A. Assume that the perturbation is small enough 
so that the system is monotone under perturbation. Monotonicity implies that, 
for all t G [0, T] and alH G f, it holds true 


4>i{t,x°,\) > 4>i{t,x°,X) > 

^i{t,x°,X) > (j)i(t,x°,X) > (l)^it,x°,X) 


so 

\4>^{t,x°,X) - 4>^{t,x°,X)\ < 4>i{t,x°,X) - 4>iit,x°,X) (52) 

In addition, [31 Lemma 1] implies that 

||0(t,x,A) - (()(t,y, A)||i < ||a; - y\\i, Vt G [0,T] (53) 
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Therefore, 


\\(l){t,x°,X) - A)||i 

< \\(l>it,x°,X) - A)||i 

< A) - A)||i + A) - x°, X)\\i 

< mt,x°,X) - + llx^'^(A) - x^^iX)\\i 

+ lk'''’(A) - A)||i + \\(j){t,x°,X) - (j3{t,x°,X)\\i 

< p° - a;"'i(A)||i + Wx^^iX) - a;®‘i(A)||i + \\x^^{X) - i°||i + p° - a;°||i 

(54) 


where the first inequality follows from (521, the second and third by triangle 
inequality, and the fourth by (53). Also, by exchanging the two terms of the 
difference after the first inequality in (54), we have 

Mt,x°,X) - A)||i < IlirO - x'^^mi + lk"‘'(A) - x'^'^(A)||i 


HA) - a 


.olli + lk°-a;°||i 


(55) 


The result follows by combining (54) and (|55|) . 


A.5 Robustness Bounds from Sensitivity Analysis 


(34) implies that 


— (ii(t) - Xi{t)) = {Xi{t) - Xi{t)) - {gi{x,a,R) - gi{x,a,R )), i€£ 

which implies that 


— {x,it) - Xi(t)) 


< |Ai(t) - Ai(t)| + \gi{x,a,R) - gi{x,a,R)\ 


Summing over i G £, and utilizing (42), we get that 


E 

i£S 


— {x,{t) - Xi{t)) 


<||A(t)-A(t)||i+L,||(i(t)-cr(t))|h 


Noting that f \xi{t) - x^{t)\ < \fAxi{t) - x^{t))\, we get 


dt 


\x{t)-x{t)\\i < ||A(t)-A(t)||i+Lg||i(t)-a;(t)||i 


Standard solution for inhomogeneous first order differential equation then gives 

(411). 


Lemma 1 For every a(t) and R{t), (42) holds true for 


Lg = 2 ("maxdHO) - min s'(a;f . 
y iGc iGc J 
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Proof 1 For every a{t), R{t), v GV, 


^ ^ I ^ ^ fij) — ^ ^ i^ij ^ ^ \di{Xi) di{xi) \ + \sj{xj) 

i^EZ 3^Et \ i^EZ 


- Xj 


< Yl d'M\xi-xi\ - Y 

ieS^ jesi' 

< Y dz{0)\xi -Xi\- Y Sjix^^Zlxj - 
i^£t 


iGEv 


Similarly, it is easy to see that, 

Y I Y diijiZ-fij) < Y d'M\xz-x^\- Y s'jZfZlZ-' 

i^EZ j^Ei' iGEZ j^Ei" 


Combining (56l and (57), we get that, for all a{t) and R{t), 
\\g{x,a,R) - g{x,a,R)\\i 

— Y1 I Y1 dijiifji — fji) + Y^ I Yh ~ fij) 

ie£ j^£- ie£ j^£ + 

= E E I E - /.o| + E E I E 

i^£+ jg£- veVi^£- jg£ + 

< 2 (maxdi{0) - mms'(a;f“)^ ||i - x||i 



I (56) 

Z (57) 

-Z) 
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